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The metallic surface state of a topological insulator (TI) is not only topologically protected, but 
exhibits an active property on curved surfaces. For an electron in the surface state of a spherical or 
a cylindrical TI (TI nanoparticle or nanowire) a pseudo-magnetic monopole or a fictitious solenoid 
is effectively induced, encoding the geometry of the system. This active property associated with 
spin-to-surface locking is mathematically formulated as spin connection, generating typically a Berry 
phase 7T. Here, we point out and demonstrate that the noninvasiveness of the surface state stems 
indeed from this Berry phase tt, the active nature of the surface state. 

PACS numbers: 73.20.-r, 73.22.-f 



Neither being a metal nor an insulator, the topolog- 
ical insulator has now been recognized as a basic form 
of solid that exhibits both gapped bulk and gapless sur- 
face states [Il-Q. Such a classification is well-defined in 
the continuum limit, while the situation is less trivial in 
the case of lattice models often employed as a concrete 
implementation of topological insulators, there remain- 
ing a question, "where actually is the surface?" A lattice 
model is sparse, and in a somewhat extreme point of 
view, existing only on sites and links, so that its surface 
is not restricted to the macroscopic boundary of the sys- 
tem, but could be also chosen e.g., such that it is partly 
extended to a rectangular-prism-shaped region (RPSR) 
as depicted Fig. [T] penetrating into the bulk. Or one can 
also think of an atomic scale closed surface isolated in 
the bulk However, in reality the protected surface 
state appears only on its macroscopic surface, exhibiting 
no symptom of penetrating into the bulk even in the case 
of sparse lattice systems. 

Why is the surface state noninvasive into the bulk? 
What prevents it from penetrating into the sparsely filled 
interior of the lattice models? In this Letter we demon- 
strate that the existence of a Berry phase 7r, or a spin 
connection associated with what is often called spin-to- 
surface locking [5Uloj|. plays a central role in this issue. 
Though existence of a protected surface state is a defin- 
ing property of the topological insulator, topological pro- 
tection does not exclude the possibility of finite-size gap 
opening. As we have demonstrated previously 1^, 111, 
Dirac electrons on the surface of a topological insulator 
exhibits an "active" property, reflecting the geometry of 
the sample. On a cylindrical surface a fictitious solenoid 
threading the cylinder is effectively induced [loj ]. while 
in the case of a sp herical system, an effective magnetic 
monopole [ll|, Il2| is induced, determining the gapped 
electronic spectrum on the surface. 

A Dirac electron on the surface of a topological insu- 
lator, especially, its spin state is susceptible of two types 
of constraints, and "locked" both in the momentum and 
real spaces. Spin-to-momentum locking is a direct conse- 
quence of the strong spin-orbit coupling in this system. 




FIG. 1: (Color online) Which is the genuine surface? 



Here, we focus on another phenomenon that manifests on 
a curved surface, often represented by the term, "spin- 
to-surface locking" . Through the bulk-surface correspon- 
dence, the entangled nature of the spin and the config- 
uration spaces encoded in the bulk Hamiltonian is tran- 
scribed to the surface Dirac equation. The helical surface 
state thus inherits a geometrical constraint imposed on 
its spin state, and an electron in this state is susceptible 
of a specific type of Berry phase, or the spin connection, 
inducing an effective monopole or a flux tube. In the 
somewhat special case of cylindrical geometry, the con- 
straint on spin manifests as spin-to-surface locking, i.e., 
the spin of the surface state is constrained onto the tan- 
gential plane of the curved surface [oHici] . Appearance of 
the spin connection in the surface Dirac equation is more 
universal, unrestricted to the case of specific geometry. 

There is an inverse effect in the "active" property of 
the surface state. Along a flux tube of strength tt (half 
unit flux quantum) piercing a TI sample a pair of gap- 
less helical modes bound to the tube is induced. These 
ID helical channels are shown to be perfectly conducting 
13], and topologically protected as well jjj, EH- In 
the presence of a surface at which the flux tube is termi- 
nated, how are these ID channels connected to the 2D 
helical surface states? In Fig. [5] we demonstrate that the 
noninvasive 2D surface state becomes gradually invasive 
into the bulk with the aid of the flux tube. When the 
total amount of the flux is not precisely tt, penetration 
of the surface state into the bulk is exponentially sup- 
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FIG. 2: (Color online) Penetration of the surface wave func- 
tion along a flux tube of strength $. 



pressed. When the flux is exactly ir, the surface state 
can penetrate into the bulk as deeply as the system's 
configuration allows it. In a sense the 7r-flux drags the 
surface state into the bulk, making it invasive. 

This Letter is intended to reveal the nature of this non- 
invasive metallic state that appears on topological insu- 
lator surfaces, especially highlighting its relation to the 
"activeness" of the state on curved surfaces. We start 
by simulating the behavior of the surface wave function 
along a flux tube. Then, as a complementary to this, 
we analytically establish the correspondence between the 
bulk and the surface descriptions. This is achieved in the 
second half of the paper, by employing a configuration 
in which the surface state can partly penetrate into the 
bulk. To ease analytic treatments the surface is designed 
to shape a smooth hyperbolic form, which may look like 
a "drain" [see Fig. [31 panel (a)]. Mathematically, this is 
the locus of a hyperbola depicted in Fig. [3] (b) when it 
revolves around the z-axis. In the limit of sharply edged 
hole (R —J- 0) this reproduces the situation described by 
the tight-binding model employed in the first part for 
numerical simulations. 

Let us briefly describe the model employed in the tight- 
binding simulation. The model is based on the following 
3D Wilson-Dirac type effective Hamiltonian in the bulk 

HG3, 



bulk 



m(p)T z + A(p x a x + p y o-y +p z a z )T x , (1) 



where m(jp) = tuq + 1TI2P 2 are Einstein and Newtonian 
mass terms encoding a band inversion due to strong 
spin-orbit coupling. Note that two types of Pauli ma- 
trices rr and t represent physically real and orbital 
spins. It is then implemented on a cubic lattice with 
nearest-neighbor hopping terms. Periodic boundary con- 
ditions are applied in the x- and y-directions (no sur- 
faces on the corresponding sides), while the model is 
restricted in the z-direction to < z < N z — 1. We 
consider a system of N x x N y x N z and introduce a 
pair of flux tubes piercing RPSRs respectively, in the 
z and — z-dircctions at (x, y) = (^f- ~ \, — and 

at (x, y) = (j-rf - — 5, — . The actual simulation is 
done in a system of size, (N x ,N y ,N z ) = (10,20,20), in 



which a moderate strength of potential disorder is also 
included 18j. Both the 2D surface and ID helical modes 



are shown to be robust against disorder. 

Depicted in Fig. [2] is the evolution of the profile of the 
lowest energy surface wave functions when a magnetic 
flux of different strength $ is introduced. As the flux 
approaches 7r, the surface state tends to penetrate into 
the bulk along a RPSR (compare different panels of Fig. 
[2] in which only a half of the system is shown) . When 
the flux is null, the RPSR is empty. Yet, one can still 
hypothesize an electronic motion bound to it. But then, 
its energy levitates because of the spin Berry phase ir; re- 
call half-odd integral quantization of the orbital angular 
momentum. Here, since the circumference of the RPSR 
is atomically small (— 4ao with a being the lattice con- 
stant), the corresponding energy scale of finite-size quan- 
tization is huge. Clearly, he is no longer compatible with 
the gapless (zero-energy) surface state. The gapless sur- 
face state, in turn, does not penetrate into the bulk along 
the RPSR. As the flux is introduced, this Berry phase ir 
is either partly or completely cancelled depending on the 
amount inserted. Then, at least a small portion of ID 
state along the flux tube starts to merge with the gapless 
surface state. From the viewpoint of the surface state, a 
portion of the wave function is dragged into the RPSR. 
As a result the wave function gets accumulated around 
the aperture. This effect should be compared with the 
asymptotic behavior of the analytic formula. 

We have seen so far through numerical simulations how 
the surface state loses its noninvasive character when a 
flux tube is inserted in an aperture of the lattice model. 
We have seen that when the strength of the flux is pre- 
cisely 7r, it becomes completely invasive. These imply, in 
turn, that the noninvasiveness of the surface state stems 
from the Berry phase 7r, which is in a sense omnipresent. 
Penetration of the surface state into any hypothetical 
apertures of the lattice is banned by the existence of this 
Berry phase n. 

To reinforce the above argument we formulate this an- 
alytically in the remainder of the paper by solving a cor- 
responding electronic state on the hyperbolic surface as 
depicted in Fig. [3j To find the surface Dirac equation 
on this curved surface it is convenient to introduce a set 
of curvilinear coordinates (£,#, </>) fl9j . defined in terms 

of the hyperbolic surface: {^J x^ + y% — aj z = R 2 ; its 
cross section in the xz-plane is shown in Fig. [3J The orig- 
inal cartesian coordinates are expressed as x — r cos cj>, 
= £ cos ( 
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is an auxiliary parameter dependent on £ and 9. The 
derivatives are represented by 
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FIG. 3: (Color online) (a) Image of the "drain", (b) Cross 
section of the hyperbolic surface on the zz-plane (only the 
x > part is shown). 



where the unit vectors e^, eg, 
dard 3D polar (spherical) coordinates [21 
sents geometrically the radius of curvature of the hyper- 
bolic curve at ro = (xq, yo, zq): 



are those of the stan- 
?7(0) 



repre- 



n(e) = v^up 



R 



1 



2 V sin 3 9 cos 3 9 



(4) 



The subsequent analyses are based on the complex ampli- 
tudes of the surface wave function at the point (£, 0, <f>), 
which is vanishingly small when £ significantly exceeds 
the penetration depth. If this is much smaller than R, 
only the range of £ <C R is physically relevant. In this 
regime we focus on hereafter apparent singularities in the 
expressions of Eq. ((3]) cause no mathematical difficulty. 

With the aid of these new coordinates we deduce the 
surface Dirac Hamiltonian on the hyperbolic surface from 
the bulk effective theory In the standard procedure 
[lol Hil l this is done by restricting the space of state vec- 
tors \ip) associated with the bulk Hamiltonian -ffbuik to a 
set of surface states, i.e., those states that are localized 
in the vicinity of the hyperbolic surface. Any surface so- 
lution lib) of -ffbuik can be written as a linear combination 



of two basis solutions, |±) 



l±», 



i.e., — t/>+|+) + t/>_|— ), where tp± are (scalar) func- 
tions of 9 and <j>. With an appropriate choice of Kx t 2 
and |±)), |±) can be made indeed (two degenerate) zero- 
energy eigenstates of -ffbuik at the "Dirac point". The 
^-dependence of the wave function is determined such 
that it vanishes on the hyperbolic surface. The spinor 
part of the wave function |±)) can be chosen as 
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Notice that here we have chosen this single-valued [in 
contrast to the standard SU(2) spinor] with respect to 
4> — » cj) + 2tt. Though this is a confusing point of this 
formulation, whether the basis is double or single val- 
ued is simply a matter of choice [l(J. The 0-dependent 



normalization constant c(0) in |±) is defined as 



c(0) 



,(6) 



in which r(rj — £) is a measure of the integral associated 
with the volume integral element r(n — £)d£d9d(f). The 
same measure appears also in the evaluation of the matrix 
elements, (±|iibuik|±) (see below). 

The surface Dirac Hamiltonian H sur f is obtained in the 
spirit of k ■ p-theory [13, EH E3]. Or, in the language 
of degenerate perturbation theory this can be regarded 
as a secular equation for the coefficients, ip±{9,<j)); they 
are solutions of the Dirac equation, H sur [if) = Eip with 
i/> = (i/'+,'0-) T - We find the coefficient matrix H sul -f by 
evaluating the matrix elements (±|i?bulk|i) as 
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Notice that in Eq. ([8J the ^-derivative in D± is shifted by 
1/2, which is nothing but the "Berry phase" of amount 
7r. Since we have chosen the spinor part of wave function 
single- valued, the orbital angular momentum L z , defined 
as i^(6, (f>) — e tLz ^Z(9), takes formally an integral value, 
L z = 0, ±1, ±2, • • • . But due to the Berry phase ir the 
physical angular momentum L z = L z + l/2 becomes half- 
odd integral 0,(2^. In Eqs. © the £ -average (/) of a 
function /(£) is defined in terms of a ^-integral similar to 
Eqs. ©, i.e., 
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The effective "Dirac theory" on the hyperbolic surface 
is prescribed by Eqs. Q, © and ©■ We now attempt 
to construct zero energy solutions of this effective model. 
To ease physical interpretation of the results it is useful 
to modify one of the coordinates by using, instead a (di- 
mensionless) angle 9, a linear coordinate £(0) having the 
dimension of length such that 
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Notice that (rj(8) — (£))d8 is a line integral element asso- 
ciated with the locus of the point r*o = (xq, yo, zo) along 
a hyperbola at fixed <f>. Thus, at a large distance r R 
on the surface (xj/-plane) from the z-axis (8 <C 7r/4), 
—((9) can be identified as the radial component r of the 
standard 2D polar coordinates (r, <fi) , while in the oppo- 
site limit of 8 3> 7r/4, ((8) can be identified as z, the 
depth into the aperture. Since d(/d8 — (r(n — £))/(r), 
the off diagonals in H SUI { [see Eq. ((TJ] becomes in the 
(£, 0)-basis, 



by the replacement, D± — > 2?± = £>± ± Af-d^ \og{r) 
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How does the wave function penetrate (or not pene- 
trate) into the aperture? What happens to the Berry 
phase 7r on the surface sufficiently away from the aper- 
ture? Answers to these questions are encoded in the ex- 
plicit form of -ffsurf- Let us focus on the zero energy so- 
lutions for comparison with the result of numerical simu- 
lations. There are two of such solutions, either with spin 
up or down, ip£2 = e li±0 Z ± (C)e±, where e+ = (1, 0) T , 
e + = (0, 1) T , which satisfy, respectively, D±ipg2o = 0- 
This can be readily solved as 



z±(0 = 



MO 



■. exp 
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(13) 



where L± = L± + 1/2 [23[. In the asymptotic limit £ — > 
oo, A^ I Ag in the exponent can be readily approximated 
as^~-^y(l-(i)2^). Deep inside the aperture, C - 

z, and also, (r) ~ a + (£) and (i) ~ (oTJ ) become 
constant, therefore Eqs. (|13p decay exponentially under 
the convergence conditions: L + > 1/2 for Z + (Q and 
L_ < —1/2 for Z_(£) [24]. In this regime, the wave 
function decays exponentially as it penetrates deeper into 
the aperture, in other words, it actually barely penetrates 
the bulk (noninvasiveness). 

How about the opposite limit, i.e., on the surface as 
( —> — oo? In this limit the profile of the wave func- 
tions can be directly compared with those of the 2D 
Dirac equation solved in terms of the Bessel functions 
J n (\E\r/A) with the use of the polar coordinates (r, <j)). 
And also, we expect that the Berry phase 7r becomes 
ineffective on the surface, which seems a priori contra- 
dictory to Eqs. ([5]) and (|T2"j) . A clue to resolve this 
discrepancy is in the normalization of the wave func- 
tion. On the 2D surface, the wave function ip2D(r,4>) 
should be normalized in terms of the surface integral 
element, rdrd(f>, while in the normalization of VKC 0) 
this measure r is not taken into account. Indeed, what 
should be interpreted as the 2D surface wave function is 
■02^(^0) = . Here, the corresponding effective 

v< r (c)> 

"2D Hamiltonian" H2D for ip2D is deduced from H SUI f 



which can be rewritten as V± = ±Agd^ ± 9< ^ e + A ( / > C±, 



-A/( in the limit of 



by noticing Ag ~ A and A^ 
C -t -00, where £+ = L+, £_ = £_ + 1. The -i— 

factor in the normalization of ip2D compensates the ef- 
fects of Berry phase ir. Thus, as expected, the Berry 
phase 7r is shown to be ineffective on the flat surface 
away from the aperture. Since in the present limit, 
A^/Ag ~ l/(r) and (r(Q) ~ the ('-integral in the 
exponent of Eqs. (|13[) diverges logarithmically, implying 
z±{C) 



oc 



|C| • These solutions are bounded only when 



C+ < for Z+(£), and L- > for Z-(Q. This implies, 
combined with the convergence conditions for the oppo- 
site asymptotics, that the zero energy solution is possi- 
ble only when L + = for Z+((), and when L_ = — 1 
for Z_((). In these two cases Z±{Cf) becomes constant, 
consistently with the fact only the zeroth order Bcsscl 
function J (\E\r/A) is compatible with the zero energy 
condition E = 0. 

Let us finally remark how the introduction of a flux 
tube piercing the aperture modifies the above argument. 
In the extreme case of $ = ir, the Aharonov-Bohm flux 
$ and the Berry phase ir (the shift of 1/2) in Eqs. © 
and (fT"2"|) cancel out each other. As a result, the bare an- 
gular momentum L± appears in the exponent of the zero 
energy solutions (|13[) . This modifies the asymptotic con- 
dition deeply inside the aperture (in the limit of C — > 00) 
to L + > and L_ < 0. In the opposite limit (on the sur- 
face away from the aperture) the two solutions behave 
asymptotically as z± ^ oc |CI ±,C± , i- e -, formally as be- 

fore, but C± now replaced with C± = L± =p 1/2. The two 
solutions are legitimate only when C + < and > 0. 
The only possible choice of L± compatible with these 
two asymptotic conditions is L+ = L- = 0. This sig- 
nifies that the wave function deeply inside the aperture 
stays constant in contrast to the previous case (expo- 
nential decay). The 7r-flux tube transforms the surface 
state invasive, penetrating the bulk to attain the oppos- 
ing surface. The asymptotic behaviors on the surface 
are modified accordingly, reproducing those of the Bcsscl 
function J_ 1 / 2 (\E\r /A) in the limit of E -» 0. 

The surface state of a topological insulator is always 
cited as being either gapless or metallic. Then, it is given 
its unique status as its existence being a defining property 
of the bulk exotic state. However, the exotic nature of 
the surface state itself was apt to be ignored. Its raison 
d'etre is in a sense considered subordinate to that of the 
exotic bulk. Here, in this Letter we have revealed that the 
surface state is by itself an "exotic state of matter" that 
stands on its own, by reclassifying it as the noninvasive 
metallic state. 
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